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^ I Abstract. 

f^ ■ In the first part, we introduce the notion of fractional statistics in the sense of Haldane. We 

illustrate it on simple models related to anyon physics and to integrable models solvable by 

the Bethe ansatz. 

In the second part, we describe the properties of the long-range interacting spin chains. We 
Q^ I describe its infinite dimensional symmetry, and we explain how the fractional statistics of 

its elementary excitations is an echo of this symmetry. 

In the third part, we review recent results on the Yangian representation theory which 
Q-(. emerged from the study of the integrable long-range interacting models. 
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1 Haldane's fractional statistics. 

1.1 Definition. 

Haldane [^ has recently introduced a notion of fractional statistics which is independent 
of the dimension of space. This notion is not based on the monodromy properties of the 
iV-particle wave functions, but on the way the number of available single-particle states 
varies when particles are added into the system. More precisely, consider a system with 
a total number of particles N = J2j Nj, with Nj the number of particles of the species j. 
Consider now adding a particle of the species i into the system without changing its size 
and the boundary conditions. Keeping fixed the positions of the A^ particles of the original 
system, the wave function of the new (A^ + l)-body system can be expanded in a basis of 
wave functions for the added particle. We denote by Di the dimension of this basis. The 
important point is that this dimension may depend on the numbers Nj of particles in the 
original system. Assuming that this dependence is linear, Haldane defines the ^^statistical 
interaction" through the relation [|T| : 

Clearly, for bosons the numbers of available single-particle states are independent of the 
numbers A^^, and Qij = 0. For fermions, the numbers of available single-particle states 
decrease by one for each particle added, and gij = 6ij. 

One of the ideas underlying the introduction of the generalized Pauli principle (|l|) is the 
fact that bosons and fermions can be considered on an equal footing as far as state counting 
is concerned. Indeed, for bosons or fermions, the number of states of A^ identical particles 
distributed among G accessible orbitals can be written in a unified way as : 

{D,j + N-1)\ 
Wbj - 



m {D,j - 1 



with Df,{N) = G for bosons, and Df{N) = G — N + 1 for fermions. The dimensions D^j 
are the numbers of accesible states for the A^*'* particle to be added. As it should be, Dj 
decreases by one unit each time a fermion is added. This is generalized to fractional statistics 
by assuming that the total number of states with {Nj} particles is : 



ii A^.! |A({A',}) - 1]! • '■ ' 

where Di{{Nj}) is obtained by integrating (|l|) : 

D^{{N,}) + ^g.,N, = G',, (3) 

with G? = Di{{0}) a constant, which is interpreted as the number of available single-particle 
states when no particle is present in the system. Namely, G^ are the bare numbers of 
single-particle states. 



1.2 Anyon-inspired examples. 

The first example ||l|] is a very naive example. Let us model anyons as charged particles 
carrying a magnetic flux. Consider now a collection of {Ni} anyons with charges qi and flux 
(pi on a disc through which goes a total magnetic flux (pB- If we assume that the charged 
particles do not interact, the number Di of accessible states for a particle of charge qi in the 
disc is : 

where 0o = h/e is the flux quantum. Since the anyons carry magnetic flux, each time 
we introduce an anyon in the system we increase the total magnetic flux by an amount 
equal to the anyon flux. Therefore, the total magnetic flux is linear in the anyon numbers : 
(pB — <Pb ~'<~ J2j(pj^j- As a consequence, the numbers of accessible states depend on the 
number of anyons present in the system, and we have : 

dDi _ qi<pj 



9ij 



ON, 



The quantity 6ij = T^{gij + gji) coincides with the Bohm-Aharonov phase obtained by moving 
an anyon of charge qi and flux (pi around another anyon of charge qj and flux (pj. 

The second example is based on a study of the elementary excitations of the fractional 
quantum Hall effect 0. So, we consider electrons moving on a plane in which there is a 
uniform transverse magnetic fleld B. We denote hy z = x + iy the complex coordinate in 



this plane, and by Ib = Jh/ireB the magnetic length. If the electrons are not interacting, 
the energy spectrum is described by the Landau levels. In an appropriate gauge, a basis of 
wave functions in the flrst Landau level is given by : ipj{z,z) = z^ exp(— |2p/2/^). These 
functions are peacked around the circle \z\'^ = j l"^. To conflne the electrons inside a disc of 
radius R, we impose that j Ib < tiR^. This is equivalent to imposing that j < N^ where A^^ 
is the number of flux quanta going through the disc, i.e. (pB = iiR^B = N^(pQ. If there are 
A^e electrons in the system, the fllling factor v is deflned as : 

A'e tl of electrons , . 

^<i> tt o/ flux quanta 

If the electrons are subject to the coulomb interaction, the degeneracy between the Lan- 
dau states is removed. The energy spectrum was numerically studied in ref. Q. It can be 
described as follows. Suppose that the fllling factor v is slightly less than the fraction l/m, 
namely: 

N^ = m{N,-l)+n (5) 

for some integer n. Then, there will be a gap in the energy spectrum in the flrst Landau 
level. The gap is of order e^ /Ib- There are a large number of states above the gap but 
only few below the gap. These low energy states are the accessible states for n excitations 
carrying each a unit flux quantum. These excitations are called quasi-holes. The number of 
accessible states for the quasi-hole excitations can be understood using Laughlin anstaz for 
the quasi-hole wave functions 0, which for a fllling factor just below 1/m as in eq.@ are 
defined by : 

m{z) = e~S^ '^^1'/'^^ \{{zi - z,r P{z) (6) 



where Zj are the coordinates of the Ng electrons, and P{z) is a polynomial, symmetric in 
these variables. These wave functions are ansatz for the states which are below the gap. So 
different choices of polynomial corresponds to different accessible states for the quasi-hole 
excitations. These polynomials are constrained by the fact that the electrons are confined 
in the disc; i.e. the electron angular momenta should be less than N^. Therefore, 

degree(P) + m{Ne - I) < N^ 

or equivalently, the degree of P is less than the number of quasi-holes. Hence, the number of 
accessible states to the quasi-holes is the number of symmetric polynomials in A^e variables 
and of degree less than n in each of the variables. This number is j^ ,^ . It agrees with 
the numerical simulations [Q. Comparing with eq.(|^), we see that, at fixed number of fiux 
quanta A*",^, the effective dimension of the quasi-hole Hilbert space is: 

D,hin) = 2 + -{N^-n) (7) 

m 

It depends linearly in the quasi-hole number. Eq.(|^ gives the quasi-hole statistical interac- 
tion : 

gqh = T^— = — 8 

on m 

These examples illustrate how low energy collective excitations may possess a statistics 
with a fractional character. 

1.3 The Bethe ansatz and fractional statistics. 

We now give another example based on integrable models f^ . We will describe how the Bethe 
ansatz equations can be reinterpreted in such way that they code a statistical interaction 
between particles, assuming that particles of different momenta belong to different species. 
Although it is useful in revealing the fractional statistics of the elementary excitations of the 
exactly solvable models, this remark does not provide a new way of solving the Bethe ansatz 
equations. 

In the Bethe ansatz approach to integrable models, the information is encoded in the 
two-body S'-matrix. We denote it by S{k) with k the relative momentum of the scattering 
particles; we have S{k) = — exp{—i6{k)), where 6{k) is the phase shift and it is odd in k. The 
eigenstates of the A^-body hamiltonian (with periodic boundary condition) are labeled by N 
(pseudo-) moment a {fc,.} (r = 1, , ■ ■ ■ , A^), which are solutions of the Bethe ansatz equations : 



e 






n S{ks - K) for all r , (9) 



s^r 



with L the length of the system. The energy E of the eigenstate | fci, ■ ■ ■ , k^) is E = J2r e°(^r), 
where e°(fc) is some universal function, e.g. e°(/c) = k'^. 

The information about the statistical interactions is encoded in the Bethe ansatz equa- 
tions (Qj, if it is rewritten in appropriate form. Indeed, taking as usual the logarithm of 
eq.(|^), the Bethe ansatz equations become : 

Lkr = J2^i^r - ks) + 2n Ir , (10) 



with 6{k) the phase shift. Here {Ir} (r = 1,- ■■ ,N) is a set of integers or half integers, 
depending on N being odd or even, which one may choose to serve as the quantum numbers 
labehng the eigenstates, instead of the momenta {kj.}. For simphcity, let us consider the 
cases when these J^'s have be chosen to be all different, (this is the usual situation). The 
ground state corresponds to an equidistribution of the integral quantum numbers {/,} in 
an interval centered around the origin : Ir+i — Ir = 1- The excited states correspond to 
particle/hole excitations in the integral lattice for the {Ir}'- 

Ir+1-Ir = l + M!; , (11) 

where M^ is the number of holes, i.e. unoccupied integer numbers, between J^+i and Ir- The 
description of the states by the quantum numbers {Ir} is a fermionic description since there 
cannot be two integers taking the same value. In this decription, the statistics is simple, but 
the energy is a complicated function of the {Ir}- 

We now change to the momentum description in the t her mo dynamical limit, N -^ oo 
at fixed density D = N/L. It is convenient to introduce the variable x = r/N which varies 
from zero to one. The pseudo-momenta kr, the integers Ir and the numbers of holes Mj. are 
all functions of x. As usual, we define the density p{k) of particles of pseudo-momentum k 
and the density of holes Ph{k) by : 



pik) 



1 



IjirCi^^l /b^ 




PH{k) = ^;il ^\\ =p{k)M{k) 

±j[rhr~[-l_ '^r 



^r I 



The density of particles is normalized by J dkp{k) = D. The energy is then given by : 

E 



dkp{k)e\k). (12) 

Thus there is no interaction energy between particles of different pseudo-momenta. Contrary 
to the fermion description, in the momentum description the dynamics looks simple but as 
we will see, the statistics is non-trivial. 

The statistical interaction is hidden in the Bethe anstaz equations. Taking the difference 
of eqs. (p!0|) for r and (r -|- 1) in the limit iV — >■ cx) gives : 

/+00 fjU' 1 

_^ —<P{k,k')p{k')=p'{k) = - (13) 



with (j){k, k') the derivative of the phase shift: 0(fc, k') = 6'{k — k'). Eq.(|T3D is the well known 
thermodynamical limit of the Bethe ansatz equations. A close comparison with the state- 
counting equation (H) shows that they are identical provided that the following identification 
are made : 

GVL ^ P°(^) = ^ (14) 

NJL ^^ p{k) (15) 

D,{{N,})/L ^-^ ph{k) . (16) 

with the discrete sum replaced by the integral over momentum. This last idenfication is 
quite natural: the hole density ph{k) clearly represents the density of available states for 



an additional particle to be added. Furthermore, from eq.(|l^) we derive the formula of the 
statistical interaction in the momentum description : 

g(k,k') = 6ik-k') + :^ <l)ik,k'). (17) 

This shows that the dynamical interaction, which is summarized in the two-body phase shift, 
is transmuted into a statistical interaction. Similarly as the boson-fermion equivalence in 
one-dimension, Eq. (|T^ provides an illustration of a phenomenon familiar in one- dimension: 
the translation of a dynamical equation into a statistical equation. 

A simple example is provided by the Calogero-Sutherland model 0. This is a model of 
particles interacting through a 1/r^ potential. It is integrable, and its two-body S'-matrix 
is S{k) = —exp[—in{\ — l)sign{k)], with A the coupling constant. Therefore, (f){k,k') = 
27r(A — l)6{k — k'), and the statistical interaction is : 

g{k,k') = X6{k-k') . (18) 

The Bethe ansatz equation then reads : 

^ = pu{k) + \p{k) . (19) 

The bare energy is e^{k) = k"^. The coupling constant A governs the statistical interaction : if 
the density of particles of momentum k increases by a unit, then the holes density decreases 
by A. Eq. (|T8|) shows that the statistical interaction is purely between particles with identical 
momenta. In this respect, the Calogero-Sutherland system appears clearly as an ideal gas 
of particles with a fractional statistics. This property of the Calogero-Sutherland model is 
also apparent in its recently computed correlation functions [Q. Finally, from eq.(^) we 
see that the duality A < — > 1/A 0], which exchanges the system with coupling A and 1/A, 
corresponds to exchanging particles and holes. 

1.4 Thermodynamics. 

Knowing how to enumerate states, it is then possible to study the thermodynamics. In the 
thermodynamic limit, the numbers of particles {Nj}, as well as the bare numbers of available 
states {G^}, become infinite. But the occupation numbers rii = [Ni/G'^) remain finite. The 
entropy is S* = kB^ogW with ks the Boltzman constant. A notion of generalized ideal gas 
was introduced in ref . |jlO| . By definition, a system is called a generalized ideal gas if (a) its 
total energy with {Nj} particles is simply given by 

E = T.N, e? (20) 

j 

with constant e°, and if (b) its states are counted according to eq.(0). 

For such gases, the thermodynamic potential Q = —PV at equilibrium can be evaluated 
by minimizing Q = E — TS — J2j ^jf^j with respect to the variation of the densities n^. 
Here T is the temperature and /Xj the chemical potential for the species i. The resulting 
thermodynamics was described in ref. 0] : 



n = -ksTY.GUogC-^^) , (21) 

Wi 



where the functions Wi are determined by the equations : 







.og(l+u,)+i:ft.log(^)^l^. (22) 

The quantities Wi possess a clear physical meaning: they are equal to the mean value of 
accessible states per particles of species i, Wi = Di{{Nj})/Ni. 

These relations completely specify the thermodynamics of the generalized ideal gas. The 
other thermodynamical quantities can be derived from the relation : 

dn = -SdT - J2 Nidfii - PdV . 

i 

In particular, the occupation numbers rii are obtained from riiG^ = —dQ/d^i. It gives : 
Hi = J2j {B^^)ijy where S is a matrix with entries : B^ = Wi6ij + hij, with C^hij = gijG^. 

The simplest example considered in |T^ is the ideal gas of particles with a diagonal 
statistical interaction; i.e. gij = gSij, fii = fi. In this case, the statistical distribution of rij is 
then given by n, = l/(wi + g), with Wi satisfying eq.(^) which now becomes: 



Wi'^ 



;i + ^,)i-^ = exp((e°-/i)ABr) . (23) 



For g = (or g = 1), we recover the bosonic (or fermionic) occupation numbers. Eq.(|23D 
possesses a. g ^^ 1/g duality : Wi{T;g) = Wi{—T/g;l/g). This is the analogue of the 
particle/hole duality of the Calogero models which we mentioned in the previous section. 

For (7 7^ 0, 1, a special example with only one energy level, is that of anyons in the lowest 
Landau level. In ref. |TT|, the anyon thermodynamical potential was computed in the strong 



magnetic field limit using a diagrammatic expansion; in ref. |jTO[, it was derived from a state 



counting approach similar to those described in section (1.2). 

The Bethe ansatz solvable models also provide examples of generalized ideal gas, since 



in the momentum description the energy is given by eq.(12), (i.e. there is no interaction 
between the particles of different momenta), and the entropy is S* = A;^ logVT with W given 
by eq.(|^) with the correspondence (|15|,|16D. The thermodynamics of such gas is called the 
thermodynamic Bethe ansatz (TBA). It was developed by Yang and Yang in ref.|^. The 
equations governing the TBA are eqs. (|2T|j2^ ) with the correspondence (|14|,|15|,|16]) . 



2 Long-Range Interacting Models. 

We now present another model whose elementary excitations obey a fractional statistics. It 



is the XXX spin chain with long range interaction introduced by Haldane and Shastry [12 
see also [0. This is a variant of the spin half Heisenberg chain, with exchange inversely 
proportional to the square distance between the spins. It possesses the remarkable properties 
that its spectrum is highly degenerate and additive, and that the elementary excitations 
are spin half objects obeying a half-fractional statistics intermediate between bosons and 
fermions. 

There is a large family of integrable long range interacting spin chains which are defined 
as follows. We consider a spin chain with A^ sites, labeled by integers z, j, ■ • • ranging from 
1 to A^. On each sites there is a spin variable (Tj which takes two values: o"i = ±. The 
hamiltonians, which are all su(2) invariant, are of the following form : 

H = Y. hi, {P., - 1) (24) 



where Pij is the operator which exchanges the spins at the sites i and j. For translation 
invariance hij = h{i — j). Demanding the integrabihty of the model selects the functions h. 
The possible choices are : 



h{x) 



(sinh7a;)2 ' 



\^' 



[V{x) 



hyperbolic model (7 real) 
trigonometric model 
elliptic model. 



where V{x) is the Weierstrass function. When 7 -^ 00, the hyperbolic model reduces to 



the Heisenberg spin chain: 



hij 



6. 



ij+'^ 



+ 5j.j+i, and for 7 — > 0, the interaction becomes the 



1/x^ exchange. The hyperbolic model has not been completly solved for general 7, although 
a partial list of eigenstates is known. The elliptic model is even more intriguing since it 
interpolates between the Heisenberg spin chain of finite length and the trigonometric model 

0- 

The Haldane-Shastry spin chain is the trigonometric model. In the thermodynamical 

limit, A^ — i> 00, it reduces to the 1/x^ exchange model, but it also possesses remarkable 

properties at finite A^. Notably, its hamiltonian commutes with an infinite dimensional 



algebra whose two first generators are |15 



Qo 



Qi 



Es^ 



2^ cotg 



N 



tJi X tj j 



(25) 



(26) 



with Si the spin operators acting on the site i. The first generators are the usual su(2) 
generators. Together with the second ones, they form a representation of the su(2) Yangian, 
(which is a deformation of the su(2) current algebra, see section 3 for an introduction to the 
Yangians). This infinite dimensional symmetry is at the origin of the large degeneracy of 
the spectrum. The fact that the hamiltonian is Yangian invariant at finite N is particular to 
the Haldane-Shastry spin chain; in the Heisenberg spin chain, the Yangian symmetry only 
appears in the thermodynamical limit. 

In order to grasp the rules describing the spectrum, we first construct few eigenstates. 
Clearly, the ferromagnetic vacuum |r2) = | + + ■■■ + +) is an eigenstate : its energy is zero. 
The eigenstates in the one-magnon sector are the plane waves :|fc) = Y^n e^p{i'2'n'kn/N)cr~\Q), 

with pseudo-moment un k, 1 < k < (N — 1): the one-magnon energy is e{k) = f-^j k(k — 
N). In the two-magnon sectors, i.e. for states of the form |^) = J2n,m4'n,mO'n ^ml^) ^ the 
eigenstates which are not degenerate with the zero or one-magnon eigenstates are labeled by 
two pseudo-momenta fci, k2, with 1 < ki, k2 < (A^ — 1). They are given by : 



^S;;'^^ = ih - h) {u^'^' 



■mk2 I ■mki+nk2 
-r UJ 



UJ'^ + UJ'' 



uj'- 



UJ'^ 



(cu"'^^ 



+mk2 



UJ 



mk\+nk2 



with UJ = exp{i2'n' / N) . Note that these wave functions vanish if ki = ^2 but also if |A;i — ^2! = 
1. The energy of l^t'^i-'^^]) is E = e{ki) + 6(^2). 

From the two-magnon computation we learn two properties of the spectrum : (i) it is 
additive, e.g. the two-magnon energy is the sum of the one-magnon energies, but (ii) the 



8 



pseudo-momenta satisfy a selection rule : they are neither equal nor they differ by a unit. 



These rules are the general rules, and the full spectrum can be described as follows [jl6|. 
To each eigenstate multiplet is associated a set of pseudo-momenta {kp} which are non- 
consecutive integers ranging from 1 to (A^ — 1). The energy of an eigenstate |{A;p}) with 
pseudo-momenta {kp} is: 



7r^2 



H\{K}) = (j:<kp)j\{kp}) with <k)=^-j k{k~N) (27) 

Furthermore, the degeneracy of the multiplet with pseudo-momenta {kp} is described by 
its su(2) representation content as follows. Encode the pseudo-momenta in a sequence of 
{N — 1) labels or 1 in which the I's indicate the positions of the pseudo-momenta; add two 
O's at both extremities of the sequence which now has length (A^ + 1). Since the pseudo- 
momenta are neither equal nor consecutive, two labels 1 cannot be adjacent. The sequence 
corresponding to the ferromagnetic vacuum is a line of 0, those of the one-magnon states 
have A^ label and only one label 1, and so on. A sequence can be decomposed into the 
product of elementary motifs, which are series of {Q + 1) consecutives O's. The multiplicity 
of the spectrum is recovered if to each elementray motif of length (Q + 1) we associate a spin 
Q/2 representation of su(2). The representation content of the full sequence is then given 
by the tensor product of its elementary motifs. 

The magnons are the excitations over the ferromagnetic vacuum; the excitations over the 
antiferromagnetic vacuum are conveniently described in terms of spinous. For A^ even, the 
antiferromagnetic vacuum corresponds to the alternating sequence of symbols 010101 ■ ■ ■ 010. 
The excitations are obtained by flipping and moving the symbols and 1. Let us give the 
sequences corresponding to the first few excitations over the antiferromagnetic vacuum, (for 
concreteness we choose A^ = 10) : 

01010101010, antiferromagnetic vacuum (o) 
10 10 1 Oa;Oa;0 10, a two-spinon excitations (2a) 

1 0^.0 1 0^.0 10 10, a two-spinon excitations (2b), etc... 

We have inserted a x between any two consecutive labels 0. These crosses represent the 
spinon excitations, their number is the spinon number. Note that there is no one-spinon 
excitation for A^ even. By convention, we will say that consecutive crosses not separated by 
any label 1 correspond to spinous in the same orbital, while crosses separated by labels 1 
correspond to spinous in different orbitals. From the rules described above, it follows that 
the degeneracy of the excitations (2a) and (2b) are different : it is three in the case (2a) 
and four in the case (2b). These degeneracy are recovered by giving a su(2) spin half to the 
spinous and by assuming that spinous in the same orbital are in a fully symmetric states. 
Hence, in the case (2a), there are two spinous in the same orbital and therefore they form a 
spin one representation of su(2), and in the case (2b), the two spinous are in two different 
orbitals and therefore they form a su(2) representation isomorphic to the tensor product of 
two spin half representations of su(2). The fact that the spinous are spin half excitations 
can also be seen by looking at the excitations of a spin chain of length A^ with A^ odd. 

This description of the states generalizes to the full spectrum. We can classify the se- 
quences by their number M of pseudo-momenta. The spinon number Ngp of a sequence is 

9 



then defined by M = ~^ "^ . Since M is an integer, (N—Nsp) is always even : this means that 
the spinons are always created by pairs. A sequence of pseudo-momenta {kp] p = 1, ■ ■ ■ , M}, 
in the Ngp spinon sector, can be decomposed into {M+1) elementary motifs where, as before 
an elementary motif is a series of consecutive 0. We call the elementary motifs the accessible 
orbitals to the spinons. At fixed Ngp, there are Norb = (H — ^-^) orbitals. Hence, a sequence 
of pseudo-momenta {kp} corresponds to the filling of the Norb orbitals with respective spinon 
occupation numbers Up = {kp+i — kp — 2), with ko = —1 and /cm+i = A^ -I- 1 by convention. 
Since an elementary motif of length {Q + 1) corresponds to a spin Q/2 representation of 
su(2), the full degeneracy of the sequences is then recovered by assuming that the spinons 
are spin half objects which behave as bosons in each orbitals. 

The spinons are not bosons but "semions" : they obey a half fractional statistics. This 
follows from the fact that the number of available orbitals varies with the total occupation 

number [Q. Indeed, at spinon number Ngp, the number of orbitals is Norb = (1 H ~2~^)' 

Therefore, we have the statistical interaction : 

dNorb -. 1^ 
The fractional statistics of the spinons is also apparent in the spin-spin correlation function 



1^ . In the following section, we will describe how the fractional statistics of the spinons is 
encoded in the Yangian representation theory. 

The spinon description of spectrum is very similar to the description of the excitations 



of the Heisenberg spin chain given by Faddeev and Takhtajan 

Note that the model is gapless. Its low energy properties belong to the same universality 
class as the Heisenberg model. The low energy, low temperature, behavior is described by the 
level one su(2) WZW conformal field theory. The spinon formulation of the Haldane-Shastry 
spin chain provides a new quasi-particle description of the states in the WZW model ||T9[| . 

3 Algebraic Solution of the Long-Range Interacting 
Models. 

In this section we review few of the new results on integrable models and on the Yangian 
representation theory which emerged from the study of the long-range interacting models. 



But we first need to recall standard result concerning the algebraic Bethe ansatz, cf e.g. pO|. 



3.1 Algebraic Bethe ansatz and Yangians. 

We introduce the basic notion of the algebraic Bethe ansatz, using the quantum Heisenberg 
chain as an example. We consider a chain of length A^: on each site there is a spin variable 



relations 



aj. We denote by S'^\ a,b = 1, 2, the spin operators satisfying the su(2) commutation 



nab Qcd 



Sjk{ S''' Sf - S'"^ Sf ) {21 



The Heisenberg hamiltonian is : 

N 

h=j:i: s^'sti, = y: iPk,k+i - 1) (29) 



N N 



k=l ab k=l 

10 



Here, we have assumed periodic boundary conditions. As is well known, in order to preserve 
the integrability the spin operators S^'' should act on the spin half representation of su(2). 
So, the spin variables take only two values, ctj = ±, and the operator S^^ which acts only 
the j^^ spin is represented by the canonical matrix |a)(&|. 

The algebraic Bethe ansatz goes in few steps. 

•) The first step consists in constructing the local monodromy matrices Tj{u). These 
matrices are 2x2 matrices whose elements T°'^{u) are operators. The matrices Tj{u) are 
defined by : 

Tf{u) = u5^^ + \Sf (30) 

where m is a complex number, called the spectral parameter, and A a coupling constant. 
Note that the matrix Tj{u) only acts on the j*^ spin. The important point is that we can 
compute the commutation relations between its matrix elements. These relations can be 
gathered into the famous relations of the algebraic Bethe anstaz, see e.g. pO| : 

R{u - v){T{u) ® 1)(1 ® T{v)) = {l® T{v)){T{u) ® l)R{u - v) (31) 

where R{u) is Yang's solution of the Yang-Baxter equation, R{u) = u — X P, with P the 
exchange operator P{x ®y) = y®x. 

••) The second step consists in constructing the complete monodromy matrix, which 
we denote by T{u). It is obtained by taking the ordered product of the local monodromy 
matrices. Namely, 

T^\u) = Y. T^'"\u)T:^^^''{u) ■■■ T^''\u) (32) 

It admits an (-)-expansion: 



a2---ajv 



U 



-N ya.(^) = 5"'' + - (y. sA + ^ I E E sfst' 1 + 



\ fc / \j<k a 



The crucial point is the fact that the complete monodromy matrix (|52|) satisfy the relations 
(pT]) if the local monodromy matrices do. These relations are equivalent to the following 
quadratic commutation relations : 

{u-v)[ T''\u) , T'\v) ] = A (t'\u)T''\v) - T'\v)T''\u)) (33) 

An important consequence of the relations (|3T|) is that the transfer matrix T{u), which is 
the trace of the monodromy matrix, T{u) = tr{T{u)) = T^^{u) + T'^'^{u), is a generating 
function of commuting hamiltonians : 

[ nu) , T{v) ] = 

The Heisenberg hamiltonian is recovered by expanding the logarithm of the trace to first 
order: H oc du log T(u) 

Another generating function of commuting quantities is given by the quantum determi- 
nant detqT('u). It is defined by : 

det,r(M) = T^\u - X)T^\u) - T^\u - \)T^''{u) (34) 
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It commutes with all the matrix elements of the monodromy matrix : det qT{u),T{v) = 0. 
The quadratic algebra (^) is called a su(2) Yangian |^. More precisely, consider a 
T-matrix satisfying the commutation relations (|3TD or (p3D, and normalized to have a quan- 
tum determinant equal to one: detgT{u) = 1. Assume that the T-matrix possesses a (-)- 
expansion as follows : 



T'^'iX) = 6 



ab 



n=0 



-n—lj^ab 



(35) 



Then, the su(2) Yangian is the associative algebra generated by the elements tf^y For these 
elements, the relations (B3) are equivalent to : 



j.ab ±cd 

''(n+1)' ''(m) 



j.ab j.cd 



±ab j.cd 

''(n)5 ''(m+1) 



^cbA.ad f:adj.cb 



(36) 



\{tfL^tZ^~tZ^t 



cb ±ad ±cb ^ad 

(m)''(n) ''(n)''(m) 



Note that with the quantum determinant constraint, the (-)-expansion of the monodromy 
matrix can be reconstructed from its two first components t?QX and t^\^. The relations (|36D 
clearly shows the Yangians are deformation of loop algebras. 

• • •) The next step consists in diagonalizing the transfer matrix. The algebraic Bethe 
ansatz provides a way to perform this diagonalization inside a finite dimensional irreducible 
representation of the su(2) Yangian. Similarly as for the unitary representations of su(2), 
any finite dimensional irreducible Yangian representation is specified by an highest weight 
vector \VL). It is characterized by the following equations : 



T{u)\n) 



fiiu) 





f2iu) 



\Q) 



(37) 



where fi{u) and /2(m) are C- number functions, not operators. The product of these functions 
is related to the quantum determinant by : detgT(-u) = /2(m — \)fi{u). Due to the fact that 
the quantum determinant commutes with the T-matrix, only the ratio fi{u)/f2{u) encodes 
the data of the representation. Moreover, the Yangian representation with highest weight 
vector \Q) is finite dimensional if and only if this ratio satisfies |^ : 



Mu) Piu + \) 



f2{u) 



P{u) 



{31 



for some polynomial P{u). These polynomials are called Drinfel'd polynomials. The condi- 
tion (^) is the analogue of the fact that finite dimensional su(2) representations correspond 
to half integer spins. 

All the states in an irreducible Yangian representation are obtained by iterative actions 
of T21(m) on |fi) : 



^)=T'\u,)T'\u,)---T'\uM)m 



(39) 



The Bethe states, which are eigenstates of the transfer matrix, are of this form, but for 
particular values of the parameters Up. The relations determining these UpS are called the 
Bethe ansatz equations. They can be summarized as follows. Let us define a polynomial 
Q{u) of degree M whose roots are the Wp's : 



Q{u) 



M 



Ur. 



(40) 
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The state (|39D is then an eigenstate of the transfer matrix T(m) if the roots Up of Q{u) are 
such that this polynomial is solution of the following difference equation : 

t{u)Q{u) = h{u)Q{u - A) + f2{u)Q{u + A) (41) 

where t(u), a polynomial of degree N , is the eigenvalue of the transfer matrix on the Bethe 
state (PP|). Notice that eq. (^Tp at the same time gives the equations determining the Bethe 
roots Up and the eigenvalue t{u). Eq.(^T]) was introduced by Baxter in its solution of the 
8- vertex model [Q. 

Following an idea due to Sklyanin p3[, the Bethe eigenstates can then be rewritten in 



terms of the polynomial Q{u). Since the operator T^^(m) is a polynomial of degree {N — 1), 
let us assume that we can factorized it as follows, 



N-l 



T^\u) = \S-. \l{u-xu) (42) 



fc=i 



where the Xk are operators and S = J2j S^ ■ It follows from the relations (|33|) that the Xk 
are commuting operators. The Bethe eigenstates (^) are then given by 

1^) = iS-f' Q{xi) Q{x2) ■ ■ ■ Q{xN^i)\9) (43) 



The eqs.(OH3) reflect the separation of the variables, since the eigenstates are determined 



from the solutions of one equation for a function of one variable only. 

Finalizing the solution of the models consists in analyzing the Bethe ansatz equations 
and their solutions. This can analytically be done only in the thermodynamical limit, along 
the lines outlined in section 1.3. 

3.2 Quantization of the spectral parameter. 

The long-range interacting models cannot be solved using the algebraic Bethe ansatz. This 
follows from the fact the hamiltonian commutes with the T-matrix, and therefore non- 
degenerate eigenstates cannot be obtained by iterative action of the lowering operators 
T^^(u). Nevertheless, the tools of the algebraic Bethe ansatz are useful for constructing 
integrable long range interacting models and for deciphering the symmetries of these mod- 
els. 

To illustrate this fact, we now consider su(2) generalizations of the Calogero- Sutherland 
models. These models describe M particles interacting by long range forces. Their positions 
are parameterized by complex numbers Zi, i = 1,---,M, and each particle carries a spin 
0" = ±. The Hamiltonian is : 



M 



ZiZrj 



Hn = U^Af - E m, + A)7-^^ (44) 

j = l i^j \^i ^j) 

where A is a coupling constant and Pij exchanges the spins of the particles i and j. Notice 
we recover the Haldane-Shastry spin chain in the static limit A = (X). 

The construction of these models relies on the definition a monodromy matrix in which 
the spectral parameter has been quantized. More precisely, let us consider the monodromy 
matrix (|32|) but in which the spectral parameters have been shifted to (u — Di), where the Di 
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are operators, commuting among themselves and with the spin operators. More precisely, 
we define a T-matrix by p^ : 



f "^(m) = ^ fp(u) ff'^^(u) ■■■ fj 



r'(«) 



(45) 



with 



a2---<i]V 



Z'^'iu) 



[U 



Di)6''^ + \Sf 



u - Di 



The operators D^ we consider are defined as follows [^ : 

Di = Zid^^ + \Y^ OijKij - A ^ OjiKi. 



(46) 



(47) 



j>i 



j<i 



and Kij the operators which exchange the particles at positions Zi and Zj : 



where 6'ii = — 

i^jjZ,- = ZiKij. They obey the defining relations of a degenerate affine Hecke algebra : 



D^ , D, 



Ki ,'4_1 , iJfc 



Ki^i+iDi — Di^iKi^ij^i 





if A; 7^ z, ? + 1 
-A 



(48) 



In the mathematics literature, the role of the affine Hecke algebra in this context was revealed 
by Cherednik |^5[. In the physics literature, operators similar but different to the Dj, were 
introduced by Polykronakos |2^. Notice that these relations imply that : 



M 



Ku//^m{u] 



^ij, ■ 



0, with Aa/ (n) = n(^'* ~ ^i 



(49) 



i=l 



I.e. /S.m{u) is symmetric by permutation of the particles. This property follows from 
ii'jj_i_i, {u — Di){u — -Di+i) = 0, which is valid for all i. 

Since the operators D^ commute, the T-matrix ( ^5]) satisfies the RTT relation (^TJ). 
However, the positions and the spin variables are totally uncoupled since the operators 
Di commute with the spin operators. In order to couple them, we define a projection vr 
which consists in replacing the permutation Kij by the permutation Pij after it has been 
moved to the right of an expression. One can view this projection as the result of acting on 
wave functions totally symmetric under simultaneous permutations of the positions and of 
the spins. In more mathematical words, this procedure consists in quotienting the algebra 
generated by the permutations Kij and Pij by the left ideal generated by {Kij —Pij). We use 
it to eliminate the permutations of the particles by replacing them with those of the spins. 

The transfer matrix T{u) defined by 



T{u) = 7r(f (w)) 



(50) 



will then satisfy the Yang-Baxter equation if we can replace the projection of the product 
(1 ® T{v)){T{u) 1) by the product of the projections. Since, Am{u) is symmetric under 
permutation, it is equivalent to check this property for T'{u) = /S.m{u)T{u). For this to be 
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true, T{u) applied on a totally symmetric wave function must still be a totally symmetric 
wave function. Equivalent ly, we must have: 



71 



(ir,,f (n)) = P,,n (f (n)) (51) 



Since the permutation groups are generated by the permutations -ft'ji+i and Pi,i+i, eg . (|51|) 
is equivalent to : vr (Ki^i+iTi{u)Ti+i{u)] = Pi,i+iTr (Ti{u)Ti+i{u)] , with Ti{u) defined in (^6|). 
This is garanteed if the commutation relations of the degenerate Hecke algebra ( ^5] ) are 
satisfied. Thus, the relations (^81) are the necessary relations for this T-matrix to satisfy the 
i?TT-relation. 

An alternative presentation of this T-matrix was obtained in ref. |Q : 



M / 1 \ 

T"^(m) = 5"^ + A ^ Sf (52) 

i,j=i \u-LJij 

where L is the matrix defined by : Lij = dijZjdz^ + (1 "~ Sij)^OijPij, with % = Zi/{zi — Zj). 
In eq. (|52D , the projection vr has been explicitely done. 

The immediate consequences of this construction are the following. Since the T-matrix 
(pOD satisfies the relation (E^ it defines a representation of the su(2) Yangian. As explained 



in the previous section, the relation ([31| ) implies that T{u) = tr{T{u)) is a generating 
function of commuting hamiltonian. However, T{u) is not Yangian invariant since it does 
not commute with T itself. A clever choice consists in choising the quantum determinant 
detgT(M) as the generating function of commuting hamiltonians. It is the projection of the 
quantum determinant of T(-u): 

det,T(«) = vr ( H^ + ^A (53) 

where Am(m) is defined in eq . (|49|) . The hamiltonian (0) is the M~^-term in ([53|) . It is there- 
fore Yangian invariant. This invar iance has recently been imbeded in a bigger algebra which 
is a deformation a W^o algebra [^]. The quantum determinant (^) has been diagonalized 



in ref. p3| by directly diagonalizing the operators Di. 



3.3 Application to the Haldane-Shastry spin chain. 

We now explain how the previous construction can be used to derive the fractional selection 
rules satisfied by the eigenstates of the Haldane-Shastry spin chain. 

As mentioned in section 2, the Haldane-Shastry spin chain is Yangian invariant. There- 
fore, there exists a T-matrix commuting with the hamiltonian (p^ and satisfying the rela- 
tions (|3TD . It was constructed in |^. It is the limit m, A ^ oo with x = u/X fixed, of the 
T-matrix (|52|). Its expression is : 



^ ' ^x-L' 



T'^''(a;) = r+^Sf(^3-^).. (54) 



with L'j^j = (1 — 6ij)9ijPij, 9ij = Zi/zij with Zij = Zi — Zj, and 5'f'' is the canonical matrix 
\a){b\ acting on the i*'* spin only. For any values of the complex numbers Zj, the transfer 
matrix (0) form a representation of the exchange algebra (^) with u changed into x and A 
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normalized to one. The trigonometric spin chain corresponds to Zj = u^ with uj a primitive 
j\j-th YQQ^ Qf ^i^Q unity. For these values of Zj, the transfer matrix (|5^ commutes with the 



hamiltonian (|2^. 

In the representation (^4]), the quantum determinant is a pure number for any values of 
the z/s given by : 



N 



detq T{x) 



iV T 



An(x) 



(55) 



with An{x) the characteristic polynomial of the N x N matrix with entries 6ij: Ai\i{x) 
det(a; — 9). For the Haldane-Shastry spin chain Zj = u^ and we have : 



^ A^ + 1 

^Nix) = Y[ix+— j) 



(56) 



Since the monodromy matrix ( p^ commutes with the hamiltonian, the long-range in- 
teracting spin chain cannot be solved using the algebraic Bethe ansatz. A way to solve it 
consists first in decomposing the spin chain Hilbert space into irreducible sub-representation 
of the Yangian, and then in computing the energy in each of these irreducible blocks. For 
uj^ , the representation (54) is completely reducible. Each irreducible sub- 



the values Zj 



representation possesses a unique highest weight vector |A) which is annihilated by T 
and which is an eigenvector of the diagonal components of the transfer matrix, as in eq. 



12/ 



X) 



In ref.[^, it was shown that the corresponding eigenvalues of T^^(x) and T^^(x) can be 
expressed in terms of two polynomials P{x) and Q{x) : 

P{x+1) 

(Kt -I- 1 W P(f, 

r(x)|A) 



Q{x + 1] 
Q{x) 



P{x) 



|A) 



(57) 



These polynomials characterize the irreducible sub-representations. The polynomials Q{x) 
and P{x) are not independent, since the quantum determinant (|5^) take the same value in 
any of the irreducible blocks. They should satisfy : 



An{x) = Fix) Q{x)Q{x 



(5^ 



Therefore, the roots of P{x) and Q{x) are among those of Ai~^{x). This implies that Q{x) 
factorized as : 



^ N + 1 

p=i ^ 



(59) 



where the {fcp} are integers bewteen 1 and [N — 1). The equation (^) then admits solutions 
if and only if the roots of Q[x) are not adjacent, or equivalently, if and only if the integers 
{fcp} are neither equal nor adjacent. These integers will be identified with the rapidities 
labeling the eigenmultiplets of the spin chain. 

This provides a purely algebraic way to recover the rapidity selection rule. It also shows 
that the fractional statistics of the spinon excitations is an echo of the Yangian symmetry. 
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